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Abstract 

We give a link to a page on the Web on which we deposited a 
set of eight huge Linear Programming (LP) problems for Intensity- 
Modulated Proton Therapy (IMPT) treatment planning. These huge 
LP problems were employed in our recent research and we were asked 
to make them public. 

1 Introduction 

In our recent research [1, 2] we employed projection methods to solve some 
huge linear feasibility problems that arise in Intensity-Modulated Proton 
Therapy (IMPT) treatment planning. These huge Linear Programming (LP) 
problems present a challenge to every LP solver or solution method and we 
have made the case in our research for tackling them with some elaborate 
projection methods. 

Since the papers [1, 2] are available on the Web we are not going into any 
further details about the research itself here. We only wish to mention that, in 
the language of [1]: "The main advantage of projection methods, which makes 
them successful in many real-world applications, is computational. They 
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have the ability to handle some huge-size problems of dimensions beyond 
which more sophisticated methods cease to be efficient or even applicable 
due to memory requirements. This is so because the building blocks of a 
projection algorithm are the projections onto the given individual sets, which 
are assumed to be easy to perform, and the algorithmic structure is either 
sequential or simultaneous, or in-between, as in the block-iterative projection 
methods or in the more recent string-averaging projection methods." 

The purpose of this Technical Report, which we do not intend to publish 
otherwise, is to announce the availability of our data sets for these problems 
along with detailed instructions where to get them and how to read them. 
We were asked to post these data sets publicly and we are pleased to do so. 

The data sets are accessible from: 
http: / /dig.cs.gc.cuny.edu/~wei/web/?page_id=221. 

Please contact chen.wei@mgh.harvard.edu for any related questions. 

2 A Set of 8 Huge LP Problems for IMPT 

The following is the information for each LP problem associated with TaskO, 
Taskl, . . . , Task7. For details about these tasks consult the papers where 
the research is reported [1, 2]. The data are binary, double precision, in the 
order of o, c, d, e, f, A as in 

min (o, x) subject to c < Ax < d and e < x < f. (1) 

The components of the vectors c, d, e, f can be infinite numbers (DBL_MAX 
in double precision in C language). The sparse matrix A is in column or- 
dered sparse matrix form as described in the MOSEK documentation, which 
is standard in storage of huge sparse matrix. For example, it can be used in 
CPLEX and other LP solvers as well. The four arrays of the matrix A are in 
the order of: ptrb,ptre, asub, aval. For a sample C code of how to read the 
data see the Appendix below. 
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#Rows 


# Columns 


#NonZeros 


Objective 


Data 


TaskO 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Taskl 


314,546 


13,734 


74,554,123 


MAX 


~900M 


Task2 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Task3 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Task4 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Task5 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Task6 


302,491 


13,734 


62,256,376 


MIN 


~750M 


Task7 


604,982 


13,734 


124,754,745 


MIN 


"1.5G 



3 Appendix 

The sample code to read the data in C language: 

char DataFileName [50] = "./O.dat"; 
FILE* file = f open (DataFileName , "rb"); 
if(file == NULL){ 

cout <<endl<< " **** fail to open file " <<DataFileName << " ****" 

<<endl ; 
return 0; }else{ 

cout <<endl<< " Reading a, c, d, e, f and ptrb , ptre , asub , aval 

of A from the binary file " <<DataFileName <<"..."; } 
fread(a, s izeof ( double ) , numCol , file); 
fread(c, s izeof ( double ) , numRow , file); 
fread(d, sizeof (double ) , numRow, file); 
fread(e, sizeof (double ) , numCol , file); 
fread(f, s izeof ( double ) , numCol, file); 
fread(ptrb, sizeof (int), numRow, file); 
fread(ptre, sizeof (int), numRow, file); 
fread(asub, sizeof (int), numNonZeros , file); 
fread(aval, sizeof ( double ) , numNonZeros, file); 
f close ( f ile ) ; 
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